Identifying phases of quantum many-body systems that are universal for quantum 

computation 
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Quantum computation can proceed solely through single-qubit measurements on an appropriate 
quantum state, such as the ground state of an interacting many-body system. We investigate a 
simple spin-lattice system based on the cluster-state model, and by using nonlocal correlation func- 
tions that quantify the fidelity of quantum gates performed between distant qubits, we demonstrate 
that it possesses a quantum (zero-temperature) phase transition between a disordered phase and an 
ordered "cluster phase" in which it is possible to perform a universal set of quantum gates. 



Measurement-based quantum computation (MBQC) is 
a fundamentally new approach to quantum computing. 
MBQC proceeds by using only local adaptive measure- 
ments on single qubits. No entangling operations are 
required; all entanglement for the computation is sup- 
plied by a fixed initial resource state on a lattice of 
qubits. The canonical example of such a resource state 
is the so-called cluster state 0, 0] . Although a handful 
of other universal resources have recently been identi- 
fied [1,1,(1,0, there currently exists very little under- 
standing of precisely which properties of quantum states 
allow for universal MBQC. For example, given a state 
that is slightly perturbed from the cluster state, it is not 
currently known how to determine if it is a universal re- 
source. New theoretical tools are required to identify the 
properties of potential resource states that allow for uni- 
versal MBQC. 

A useful perspective to approach this problem is to 
view the resource state for MBQC as the ground state of 
a strongly-coupled quantum many-body system. With 
this perspective, we propose that the ability to perform 
MBQC is a type of quantum order - one which can be 
identified using appropriate correlation functions as or- 
der parameters. We show that a natural choice for such 
correlation functions are the expectation values of non- 
local strings of operators that can be identified with mea- 
surement sequences for performing quantum logic gates 
within MBQC. One way of understanding MBQC is that, 
by means of a set of local measurements, it is possible to 
prepare the resource states required for gate teleporta- 
tion 0, 0, S] between distant components of the many- 
body system. The performance of the MBQC scheme 
can be characterized by calculating the fidelity of the pre- 
pared resource state with the ideal one 0. This fidelity 
will depend on a set of non-local correlation functions as a 
result of the many local measurements that are required 
to prepare the resource state. (Because the fidelity of 
the identity gate is quantified by the ability to prepare 
an entangled state between two distant qubits using local 
measurements, it is closely related to the much-studied 
property of localizable entanglement (lo|.) We show that, 



for the cluster state implementation of MBQC, the spe- 
cific correlation functions corresponding to any gate can 
be calculated, and we investigate a specific model where 
the fidelities of a gate set indeed serve as order parame- 
ters identifying a cluster phase. This result suggests the 
existence of spin systems that possess a phase for which 
any state is a universal resource for MBQC. These meth- 
ods provide new tools for identifying properties of quan- 
tum many-body systems that are required for MBQC. 

Consider the following model system. The cluster state 
on a lattice C is defined as the unique +1 cigenstate of a 
set of stabilizer operators = X lt \\ v ^ j Z v , where 
(Zfj,) is the Pauli X (Z) operator at site ijl and where 
v ~ /i denotes that v is connected to fi by a bond in 
the lattice C. The Hamiltonian H = — Ylu^c has the 
cluster state as its unique ground state ll| . Although the 
terms in this Hamiltonian are many-body interactions, it 
can be realized as the effective low-energy theory of a 
Hamiltonian consisting only of two-body terms [121 ]. 

As a model system to consider how robust is this 
Hamiltonian in the presence of local perturbations, we 
supplement it with a local field term, 



H(B) = -^ eC (K^ + BX u ). 



(1) 



representing a local transverse field with magnitude 
B. We refer to a lattice with this Hamiltonian as 
the transverse-field cluster model (TFCM), and we will 
demonstrate the existence of a single zero-temperature 
phase transition in the ground state of such models 
on both a 1-D line and a 2-D square lattice, separat- 
ing a disordered phase from a "cluster phase". Rather 
than solving these models explicitly, we explore duality 
transformations that relate these models to others with 
well-understood phases and order parameters. We then 
demonstrate that the order parameters of these mod- 
els, mapped back to the TFCM, are precisely equiva- 
lent to the correlation functions in the cluster state that 
quantify the fidelity of the identity gate (i.e., teleporta- 
tion) in MBQC. That is, the ability to perform the iden- 
tity gate over a long range serves as an order parameter 
for this phase; similar results hold for other single-qubit 
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gates as well. In addition, in two dimensions, we per- 
form a similar analysis of the two-qubit CSIGN gate, 
exp(z7r|l)(l| <X> |1)(1|), which together with our singlc- 
qubit gates yields a universal gate set for MBQC. (In 
contrast, the case with a local longitudinal field instead 
of a transverse one was investigated in [l]|; this model 
demonstrates no such phase but nevertheless can still al- 
low for MBQC for some range of parameters.) 

General properties of the transverse field cluster model. 
We first present some general properties of the TFCM 
that are valid in any dimension and on many lattices, 
before investigating one- and two-dimensional models in 
detail. An immediate observation is that this model is 
self-dual. The canonical transformation of Pauli opera- 
tors given by applying the CSIGN operation between all 
neighbouring pairs of qubits takes <-» X^, and thus 
the Hamiltonian (Q} transforms as H{B) — » BH(1/B). 
This self-duality ensures that, if this model has a single 
quantum phase transition in the range B > 0, then it 
must occur at B = 1. 

Also, consider lattices which are bipartite, meaning we 
can divide the sites into two subsets £ r and Cb, labeled 
red and blue, such that the neighbours of any site are all 
of the other colour. With this colouring, the Hamiltonian 
([1]) can be written as the sum of two commuting terms, 
H = H T + i?b , where 

Hr = -E„ e£b ^ - BE M&C ** . (2) 

with H\j consisting of the remaining terms. In the fol- 
lowing, we present mappings of _ff r (equivalcntly, H\>) in 
one and two dimensions to known models, which allows 
us to identify the phases and relevant order parameters. 

One dimension. Consider the TFCM on a 1-D lattice 
with fixed boundary conditions - a line. A state of a 1-D 
lattice cannot serve as a universal resource for MBQC; 
however, it will be illustrative to consider this model as 
a prelude for studying higher dimensions. The Hamilto- 
nian JI| on a line with boundary terms is 

H(B) = -E^ 1 (Z. t _ 1 A 4 Z 4+1 + BXi) 

— X1Z2 — BX\ — Zn-iXn — BXn ■ (3) 

The ground state of this Hamiltonian is non-degenerate, 
and for B = is given by the 1-D cluster state on a 
line. Pachos and Plenio [lj| have shown explicitly that 
this model (with periodic boundary conditions) exhibits a 
quantum phase transition at \B\ = 1, and that the local- 
izablc entanglement length remains infinite for all values 
\B\ < 1. Their method makes use of the Jordan- Wigncr 
transformation to yield a linear fermionic system. We 
provide a more direct transformation to a known model 

the transverse-field Ising model [loT ] - which provides a 
natural generalization to higher-dimensional lattices. 

Our duality transformation is as follows. On red (even) 
sites, the Pauli operators transform as 

X 2 j -> X 2j , Z 2 j -> (Y[ 3 k=1 X 2 k~i) Z 2 j ■ (4) 




FIG. 1: (a) The duality transformation of Eqs. (|4I5|I on a 1-D 
line, (b) A generalization of this duality transformation to a 
2-D square lattice. 

On blue (odd) sites, the Pauli operators transform as 

X23-1 -> X 2j -i , Z 2j -i -> Z 2 j-i(J\^ =j X 2 k) ■ (5) 

This mapping is canonical, meaning the new Pauli matri- 
ces Xj and Zj satisfy the correct commutation and anti- 
commutation relations. An illustration of this transfor- 
mation is presented in Fig. [Ha). We emphasize that this 
duality transformation is non-local, and thus the proper- 
ties of a system for MBQC arc not preserved under this 
mapping. However, as we now demonstrate, the phases 
and order parameters of this dual model are well-studied 
and will allow us to completely classify the phases as well 
as calculate the fidelities of the MBQC quantum gates in 
the original TFCM. 

We consider only the case where N is even. In terms 
of transformed Pauli operators, the Hamiltonian H T acts 
only on red sites and has the form 

H T = -Z 2 - EfJ 2 (Z 2 {i-i)Z 2i + BX 2l ) , (6) 

The Hamiltonian Hy> is similar, acting only on blue sites, 
with a Z boundary term at j = N. This mapping on the 
TFCM, then, yields two identical transverse-field Ising 
models, one on each of C r and £t>- Each has a local 
Z field term which breaks the symmetry in the ordered 
(|£>| < 1) phase and specifies a unique ground state. The 
ground state of the total lattice is then non-degenerate 
and is given by the product state of these two unique 
ground states. 
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The solution to this known model allows us, via the 
duality transformation, to completely characterise the 
TFCM. For example, the phases of the TFCM are spec- 
ified by the well-studied phases of the transverse-field 
Ising model; in particular, there is a unique quantum 
phase transition at \B\ = 1 [lq |. Also, the well-known 
order parameters for the transverse-field Ising model can 
be mapped, using the duality transformation, to order 
parameters for the TFCM. In the ordered phase of the 
transverse-field Ising model, the correlation functions 
(ZZ) (for both colors) are long ranged. (Specifically, 
Ymxk^ZiZi+k) = (1 - \B\ 2 ) 1/4 for \B\ < 1 0.) We 
can use this result to make a corresponding statement 
about correlation functions for the TFCM. By reversing 
the duality transformation, we have 



(a) 



Z 2i -iZ 2j -i) - (Zai-xdlQ-yaO^-i) = (iF^a) . 

(7) 

(8) 



(z 2i z 2j ) - (z 2t (Ui-)x 2k+1 )z 2j ) = (Ut]K 2k+1 ) 



That is, in the phase \B\ < 1 wherein (ZZ) is long- 
ranged, the string-like operators corresponding to the 
product of even (or odd) stabilizers Ki in the TFCM are 
also long-ranged, with the limiting value (1 — li?! 2 ) 1 / 4 . 
These two correlation functions are all that is needed to 
calculate the fidelity of the resource state for the identity 
gate with the ideal maximally- entangled state, and it is 
found to be > 1/4 for all \B\ < 1. (The average fidelity 
of a randomly chosen state yields 1/4.) The same calcu- 
lation for other single-qubit Clifford gates [l7| and for an 
arbitrary Z-rotation U z {9) = exp(—i8Z) (a non-Clifford 
gate), yields the same result 

Thus, this duality transformation has allowed us to 
prove our desired results: First, that the TFCM does in- 
deed possess a phase, given by \B\ < 1, which we denote 
the cluster phase. The order parameters of this phase, 
given by products of even or odd stabilizer operators Ki, 
demonstrate that quantum gates can be performed with 
high fidelity (relative to a randomly-chosen state) using 
any state within this phase. The ground states in this 
phase are indeed "robust" against variations in the pre- 
cise value of B. However, the one-dimensional cluster 
state is not a universal resource for MBQC, and so we 
direct our attention to a two-dimensional model. 

Two dimensions. We consider a square lattice; the 
cluster state on this lattice is a universal resource for 
MBQC. This lattice is bipartite, and thus we can define 
the commuting Hamiltonians H t and as above. We 
use a natural generalization of the 1-D duality transfor- 
mation, as follows. On red sites, Pauli operators trans- 
form as Xfj. -» X^ and Z^ — > (JJ^^^X^Z^, whereas on 
blue sites, X^ -> X^ and Z^ -> Z^fJl^^X^). Here, 
fj,' > (J, (p, 1 < n) denotes that fx' lies in the upper (lower) 
cone relative to /i as in Fig. HJb). Again, one can easily 
verify that this transformation is canonical. 
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FIG. 2: (a) A measurement pattern on the cluster state that 
localizes entanglement between sites ai n and a ou t, where X 
(Z) denotes a measurement in the X-basis (Z-basis). The 
two string-like stabilizers, centred on sites connected by the 
shaded red and shaded blue diagonal lines, have long-ranged 
expectation values in the |B| < 1 phase; these correlation 
functions directly quantify the fidelities of single-qubit gates 
between ai n and a out in MBQC. (b) The measurement se- 
quence corresponding to the CSIGN gate between a and b. 
The expectation of four stabilizers characterises the CSIGN 
gate: K ain K s K aaut , K hin K 4 K hollt , K X K± and K 2 K 3 . These 
stabilizers can be appended with diagonal strings of red (blue) 
stabilizers in the direction of the arrows (and terminated with 
Z measurements as in (a)) to reach distant qubits. With X 
measurements on qubits 1-4, the resulting state provides the 
CSIGN transformation. 



Under this mapping, each stabilizer maps to a 
monochromatic operator consisting only of Z terms. 
Non-boundary stabilizers map to products of four Z op- 
erators on the corners of a fundamental plaquette □; 
boundary conditions can be chosen such that boundary 
stabilizers map to two-Z and oub-Z terms. The Hamil- 
tonian H r (H^) on £ r (£h) maps to 



H 



zZz z~ B z2*v 



(9) 



plus boundary terms (not shown) which ensure a non- 
degenerate ground state for all B. This model possesses 
a phase transition at \B\ = 1 0,111. Thus, through 
this duality map, we know that the 2-D TFCM has a 
phase transition at \B\ = 1, and we use the term clus- 
ter phase to denote the \B\ < 1 phase. In addition, this 
model of Eq. ([9]) is dual to the anisotropic quantum or- 
bital compass model (AQOCM) 0, ED, H3, H3 , with a 
mapping that also locally maps the boundary terms). 
The key advantage of the AQOCM is that it contains 
only two-body terms in the Hamiltonian, and is therefore 
very amenable to numerical investigation. For example, 
the projected entangled-pair state algorithm applied to 
this model provides very strong evidence that the phase 
transition is first order j24[. The model also possesses 
correlation functions for an Ising order parameter that 
simulations indicate are long-ranged for \B\ < 1 



24] 



Inverting this duality transformation and returning to 
the TFCM, these Ising- type correlation functions map 
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onto strings of monochromatic stabilizers along diagonal 
lines in the square lattice (see Fig.[2ja)). Again using the 
correlation functions for single-qubit gates given in [9( , we 
find that these strings of monochromatic stabilizers char- 
acterize the fidelities of the identity gate and a generating 
set of single-qubit gates between two distant points, and 
serve as order parameters for the cluster phase. 

In addition, in this 2-D model we can consider two- 
qubit gates. We make use of the elementary measurement 
pattern for a CSIGN gate on two qubits which are sub- 
sequently swapped, as given in Refs. 0, Q and shown in 
Fig. [2b) . The desired long-ranged correlation functions 
on the AQOCM are of the form of 4-body correlators 
(i' (ij0 )Z(j Jii )i'( i+ljji )Z( i+ljl )), where is an intermedi- 
ate column between j and j\ . Such 4-body correlation 
functions should be possible to numerically evaluate in 
the AQOCM using recent techniques. The CSIGN to- 
gether with the above single-qubit gates yields a universal 
gate set, and thus the cluster phase is indeed character- 
ized by the fidelities of a universal gate set for MBQC. 

Discussion. Using the TFCM as an example, we have 
demonstrated the utility of correlation functions corre- 
sponding to quantum gates as order parameters to iden- 
tify a phase according to its usefulness for MBQC. The 
perspective of quantum-computational universality of a 
state as a new type of quantum order may assist in identi- 
fying new quantum systems that can be used for MBQC. 

The behaviour of the TFCM contrasts with the model 
considered in which is the cluster state Hamilto- 

nian perturbed by a local Z-ficld. In that model the 
gate correlation functions discussed in this paper become 
short-ranged at any non-zero perturbation. However, by 
pre-processing with certain local filtering operations it is 
still possible to perform MBQC for sufficiently low field 
and sufficientl y lo w temperature [l3|. Unlike the TFCM, 
the model of [13j does not undergo a phase transition. 
These behaviours are very reminiscent of the quantum 
Ising model in one dimension; where there is a broken 
symmetry that disappears at a phase transition for suf- 
ficiently large transverse field but longitudinal fields de- 
stroy the ground state order without any phase transi- 
tion. 

One could also ask whether these ordered phases per- 
sist to finite temperature. As our model is gapped ex- 
cept at the phase transition, it is possible for a finite- 
sized thermal system to be cooled to have arbitrarily high 
overlap with the ground state (although this becomes a 
challenge close to the phase transition). In one dimen- 
sion, the fact that the transverse-field Ising model does 
not maintain an ordered phase at any finite temperature 
demonstrates that the 1-D TFCM does not cither. In two 



dimensions, it is less clear. For this reason it would be 
worth investigatin g th e TFCM on a three-dimensional 
lattice such as in [ll|, for which the B = model is 
known to allow for fault-tolerant MBQC at finite tem- 
perature HE EE Sij]. 
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